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deformation  of  wet  snow  are  given  for  a variety  of  conditions.  These  results  are  in  agreement  with  existing  experi- 
mental evidence.  The  model  accurately  simulates  the  transient  nature  of  the  deformation  and  the  effect  of  water  con- 
tent on  the  quasi-stable  density  of  wet  snow  subjected  to  a constant  stress.  ^ 


PREFACE 

This  report  was  prepared  by  Dr.  Samuel  C.  Col  beck,  Geophysicist,  of  the  Snow 
and  Ice  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineer- 
ing Laboratory.  The  research  was  funded  by  DA  Project  4A161 102AT24,  Research 
in  Snow,  Ice,  and  Frozen  Ground,  Scientific  Area  03,  Research  in  Terrain  and  Cli- 
matic Constants,  Scientific  Effort  01 , Cold  Environment  Factors,  Work  Unit  001 , 
Research  in  Snow  Mechanics , and  by  DA  Project  1 T 1 6 1 1 02B52 A,  Mobility  and  En- 
vironmental Research , Task  02,  Military  Aspects  of  Cold  Regions  Research , Work 
Unit  003,  Research  in  Snow  Mechanics  (QCR  CARDS  111a,  110b). 

Dr.  Malcolm  Mellor,  Dr.  William  Harrison,  and  Nancy  Parssinen  of  CRREL  tech- 
nically reviewed  the  manuscript. 


LIST  OF  SYMBOLS 


d thickness  of  water  film 

D depth  (or  displacement)  of  melted  cap 

F force  between  particles 

k thermal  conductivity  of  ice 

L latent  heat  of  fusion 

n%  moles  of  solute 

moles  of  water 

p pressure  in  the  liquid  phase 

Po  pressure  in  the  liquid  at  the  center  of  the  contact 

Fb  pressure  at  the  bond 

Pc  capillary  pressure 

q water  flux  in  film 

Q heat  flux 

r polar  coordinate 

rt  radius  of  air  bubble 

rb  radius  of  contact 

ro  initial  particle  radius 

/p  particle  radius 

R gas  constant  for  water  vapor 

Tb  temperature  of  bond 

T0  273  K 

Tt  temperature  of  curved  particle  surface 
p viscosity  of  water 

p density  of  snow  sample  considering  ice  phase  only 
ps  liquid  water  density 

p0  initial  ice  density  of  snow  sample 

P,  density  of  ice 

o8g  liquid/vapor  interfacial  tension 
oH  solid/gas  interfacial  tension 
solid/liquid  interfacial  tension 


j! 

■ji 

SUMMARY 


The  rheology  of  wet  snow  controls  over-snow  movement  in  the  cold  regions 
during  much  of  the  year.  The  deformation  of  wet  snow  is  equally  important  to 
traffic  movement  over  highways  and  to  the  release  of  wet  snow  avalanches.  The 
thermodynamics  of  wet  snow  controls  its  rapid  deformation  in  a manner  unlike 
the  deformation  of  dry  snow.  Because  of  the  presence  of  liquid  water,  wet  snow 
can  deform  rapidly  by  melting  at  the  inter-particle  contacts  and  refreezing  at  the 
curved  ice  surfaces.  The  rate  of  deformation  is  controlled  by  the  parameters  which 
determine  the  temperature  difference  between  the  inter-particle  contacts  and  the 
curved  ice  surfaces.  The  temperature  of  the  contact  is  determined  strictly  by  the 
force  between  the  particles  and  the  size  of  the  contact.  The  temperature  of  the 
curved  surface  is  determined  by  the  particle  size,  ionic  impurity  content  of  the 
liquid  phase,  and  air  bubble  size  or  capillary  pressure.  As  the  inter-particle  bond 
size  increases,  the  inter-particle  stress  decreases  and  the  bond  temperature  increases 
until  the  temperature  differences  around  a particle  are  equalized.  At  that  time  the 
rapid  thermodynamic  deformation  described  here  ceases,  and  further  deformation 
must  proceed  by  more  efficient  packing  arrangements  and/or  the  normal  processes 
of  crystalline  deformation. 

Using  a simple  geometrical  packing  arrangement  (see  cover  of  this  report)  and 
the  principles  of  thermodynamics,  a simple  model  of  wet  snow  deformation  is 
constructed  here.  The  results  of  this  model  are  qualitatively  similar  to  the  experi- 
mental results,  including  a lower  value  of  maximum  density  at  larger  values  of 
capillary  pressures.  The  effect  of  ionic  impurities  on  the  maximum  bond  size  is 
also  calculated  and  found  to  be  the  same  as  increasing  the  capillary  pressure.  The 
hardening  effects  of  increasing  capillary  pressure  and/or  ionic  impurity  content 
have  practical  implications.  For  example,  the  optimum  use  of  salt  for  snow  re- 
moval on  highways  and  for  hardening  wet  snow  for  increased  over-snow  mobil- 
ity can  be  achieved  only  through  the  understanding  of  wet  snow  rheology. 
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THERMODYNAMIC  DEFORMATION  OF  WET  SNOW 

Samuel  C.  Colbeck 


INTRODUCTION 

The  rheology  of  snow  at  subfreezing  temperatures 
has  been  extensively  researched  and  much  useful  infor- 
mation has  been  generated  (see  Mellor  1975).  Unfor- 
tunately, the  rheology  of  snow  which  contains  liquid 
water  has  been  almost  totally  ignored.  This  disparity 
does  not  arise  because  wet  snow  is  any  less  important 
than  snow  at  lower  temperatures  — witness  the  des- 
tructive powers  of  wet  snow  slides  and  the  highway 
hazards  caused  by  the  compaction  of  wet  snow  into 
ice.  It  probably  arises  from  the  greater  difficulties 
of  handling  wet  snow  samples  and  of  understanding 
the  complex  interaction  between  the  thermodynamic 
and  mechanical  aspects  of  wet  snow  deformation. 

For  a given  stress,  the  deformation  of  wet  snow 
generally  occurs  much  more  rapidly  than  can  be  ex- 
plained by  intracrystalline  straining  and/or  more  ef- 
ficient grain  packing;  hence,  some  thermodynamic 
explanation  must  be  invoked.  Colbeck  (1973)  des- 
cribed some  of  the  important  aspects  of  the  phase 
thermodynamics  of  wet  snow.  The  equilibrium  temp- 
erature among  the  three  phases  of  water  in  wet  snow 
is  controlled  by  the  size  of  the  spherical  grains,  water 
content,  grain  contact  area,  impurity  content,  por- 
osity, and  stress.  The  thermodynamic  and  geometric 
principles  which  relate  the  deformation  to  these  vari- 
ables are  used  here  to  describe  the  densification  of 
wet  snow  at  high  strain  rates.  The  slow  deformation 
of  wet  snow  in  situations  where  intracrystalline  creep 
or  grain  repacking  are  the  dominant  mechanisms  is 
not  considered  here. 


THERMODYNAMICS 

The  thermodynamic  deformation  of  wet  snow  is 
controlled  by  the  temperature  difference  between 
the  inter-particle  contacts  and  the  ice/fluid  interfaces 


Figure  1.  Schematic  rep- 
resentation of  two  par- 
ticles moving  together  by 
melting  at  their  mutual 
contact.  The  size  of  the 
particle  ip,  size  of  the 
bond  rfc>  and  thickness 
of  the  water  film  d are 
shown. 


(Fig.  1).  The  temperature  at  a contact  is  determined 
exclusively  by  its  area  and  the  average  normal  force 
acting  across  it.  During  rapid  deformation  no  shear 
forces  can  be  supported  by  an  inter-particle  contact 
because  the  ice  surfaces  are  separated  by  a water  film 
of  sufficient  thickness  to  prevent  ice-to-ice  contact. 
The  thickness  of  this  film,  as  determined  by  the  rate 
of  deformation,  is  described  below. 

The  temperature  of  the  free  surface  away  from  tire 
inter-particle  contact  is  determined  by  the  radius  of 
the  particles,  radius  of  the  air/liquid  interface,  and  im- 
purity content  of  the  liquid  phase  in  tire  pores.  De- 
pending on  the  direction  of  heat  flow,  there  are  three 
possible  mechanical  states. 

1.  When  the  surface  temperature  7^  is  less 
than  the  contact  temperature  Tb,  heat  flows 
away  from  the  contact  and  produces  freezing 
in  the  liquid  layer.  Thus  the  inter-particle 
strength  tends  to  increase.  This  regime  oc- 
curs when  the  snow  is  frozen  either  by  spread- 
ing salt  on  a wet  snow  surface  or  by  withdraw- 
ing the  liquid.  Accordingly,  the  relatively  rapid 
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thermodynamic  deformation  ceases  and  the  wet 
snow  deforms  by  the  same  mechanisms  as  dry 
snow.  If  a higher  stress  or  more  liquid  is  intro- 
duced, the  contacts  may  melt  again. 

2.  When  the  pore  surface  and  contact  temp- 
eratures are  equal,  no  heat  flow  occurs,  but  a 
significant  inter-particle  strength  can  develop 
slowly  due  to  ice-to-ice  adhesion.  The  liquid 
tension  may  cause  some  additional  inter-par- 
ticle strength  (Moore  and  Millar  1971,  Skelton 
1975). 

3.  When  the  lowest  temperature  occurs  at 
the  inter-particle  contact  because  of  the  depres- 
sion of  the  phase  equilibrium  temperature  by 
the  contact  stress,  heat  flows  toward  the  con- 
tact causing  a continuous  melting  of  ice  at  the 
contact.  This  melting  at  the  high  pressure  con- 
tact is  accompanied  by  water  flow  away  from 
the  contact  toward  the  low  pressure  pore 
water.  Thus  water  flow  and  heat  flow  are  bal- 
anced, resulting  in  melting  at  the  contacts,  re- 
freezing on  the  curved  surfaces,  and  densifica- 
tion  as  the  particle  contacts  move  closer  togeth- 
er and  the  porosity  is  reduced  by  refreezing. 


When  the  lowest  temperature  occurs  at  the  bonds, 
rapid  densification  can  occur  even  at  moderate  levels 
of  normal  stress  transmitted  across  the  contacts.  The 
contact  temperature  during  rapid  densification  is  in- 
dependent of  impurities  in  the  pore  water  because 
meltwater  moves  away  from  the  contact.  Accord- 
ingly, the  equilibrium  temperature  7J,  for  the  nearly 
flat  contact  surface  is  given  by* 


(i) 


where  Pb  = pressure  of  the  bond 
T0  = 273  K 

L = latent  heat  of  fusion 
pa  = liquid  water  density 
ps  = density  of  ice. 

The  phase-equilibrium  temperature  on  the  rounded 
particle  surface  depends  upon  the  saturation  regime  of 
die  wet  snow  (Colbeck  1973).  The  equilibrium 


* The  effect  of  dissolved  air  has  been  neglected  because 
of  uncertainty  about  the  air  content  of  the  meltwater 
at  the  contacts. 
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temperature  in  the  funicular  regime  of  high  liquid  sat- 
uration is  determined  by  the  ice  particle  radius  rp  and 
air  bubble  radius  ra  according  to 


T = ^ T°/  1 _ 1 \qgg  _ 2T0  °sB 
L Ip,,  pj  ra  Lps  rp 

RTq  nt 


(2) 


where  oag  = liquid/vapor  interfacial  tension 
asE  = solid/liquid  interfacial  tension 
R = gas  content  in  water  vapor 
njnw  = temperature  depression  caused  by  the  pre- 
sence of  ionic  impurities. 


Similarly,  the  equilibrium  in  the  pendular  regime  of  low 
saturations  is  determined  by  the  particle  size,  ionic  im- 
purity content,  and  capillary  pressure  Pc  according  to 


T - 

s ' PaL  c‘  Lp. 


gsg  RTo  "s 

. rp  1 nw 


(3) 


where  asg  is  the  solid/gas  interfacial  tension. 

Ionic  impurities,  capillary  pressures  and  small  par- 
ticle sizes  reduce  the  temperature  of  the  curved  parti- 
cle surface;  the  heat  flow  to  the  bonds  is  reduced  ac- 
cordingly. The  rate  of  deformation,  ultimate  bond 
size,  and  maximum  density  achievable  by  rapid  ther- 
modynamic densification  are  also  reduced.  In  snow 
which  is  unstressed  but  saturated,  particle  size  in- 
creases rapidly  to  the  quasi-stable  diameter  of  1 to  3 
mm  (Wakahama  1968,  Colbeck  1973).  At  this  time 
the  effect  of  particle  size  on  the  phase  equilibrium 
temperature  is  only  10' 5 to  10*4  °C,  or  about  the 
same  as  the  effect  of  impurities  in  natural  snow. 

Capillary  pressure  increases  as  water  saturation  de- 
creases, as  shown  experimentally  by  Colbeck  (1973). 
Thus,  at  high  water  saturations  where  Pc  is  small,  the 
curved  surface  temperature  is  high;  therefore,  heat 
flow  into  the  inter-particle  contact  and  rapid  densifi- 
cation are  more  likely.  As  water  saturation  is  reduced 
and  capillary  pressure  increases,  the  rate  of  densifi- 
cation is  reduced  accordingly.  The  experimental  re- 
sults of  Tusima (1973) and  Kinosita  (1963)demon- 
strate  this  fact.  In  principle  the  effect  of  an  increase 
in  capillary  pressure  could  be  as  large  as  -0.01°C  and 
profoundly  influence  the  response  of  wet  snow  to 
stress.  In  practice,  however,  this  effect  is  greatly  re- 
duced by  the  time  required  to  achieve  large  capillary 
pressures  in  large  grained  materials. 
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Harris  and  Morrow  (1964)  showed  that  reduction 
of  liquid  pressure  is  delayed  during  drainage  because 
of  the  slow  process  of  liquid  movement  through  thin 
films  on  particle  surfaces.  Their  results  suggest  that 
Colbeck’s  (1973)  capillary  pressure-water  saturation 
curve  does  not  quite  represent  true  capillary  equili- 
brium, because  the  experiment  was  completed  too 
quickly  to  allow  complete  drainage  at  each  pressure 
level.  Thus,  while  capillary  pressure  can  have  a large 
effect  on  the  densification  of  wet  snow,  many  days 
may  be  required  for  complete  liquid  drainage  from 
wet  snow  subjected  to  large  capillary  pressures.  Much 
experimental  work  remains  to  be  done  on  the  capil- 
lary equilibrium  in  the  pendular  regime  of  saturation 
to  determine  the  circumstances  necessary  for  achiev- 
ing large  capillary  pressures  in  wet  snow  samples.  The 
main  experimental  problem  is  the  loss  of  liquid  con- 
tact between  the  sample  and  the  porous  plate  through 
which  the  liquid  is  extracted.  Occasionally  it  is  pos- 
sible to  measure  large  capillary  pressures  in  wet  snow 
(Colbeck  1976),  suggesting  that  loss  of  contact  can- 
not necessarily  be  predicted  in  advance.  This  prob- 
lem, which  is  common  for  porous  materials  with  large 
grains,  does  not  preclude  the  application  of  this  theory 
of  thermodynamic  densification  to  samples  subjected 
to  large  capillary  pressures.  However,  it  must  be  recog- 
nized that  large  capillary  pressures  are  difficult  to  ob- 
tain experimentally. 


HEAT  AND  MASS  FLOW 

When  a lower  temperature  occurs  at  the  inter-par- 
ticle contact  than  at  the  particle  surface,  the  heat  flux 
Q approaching  the  contact  can  be  determined  by 

T - T 

Q=k-±-±-  (4) 

'’b 

where  k is  the  thermal  conductivity  of  ice  and  rb  the 
radius  of  the  inter-particle  bond.  The  contact  melts 
at  a rate  given  by 

dD_  k TrTb  . . 

dt~Lp,  rb 

where  D is  the  depth  of  the  melted  cap  (see  Fig.  2). 
Equation  5 gives  the  rate  of  approach  of  particle  cen- 
ters, but  the  particle  size,  impurity  level,  water  con- 


Figure 2.  Cross-sectional  view  through 
the  center  of  a sphere,  showing  the  ini- 
tial radius  r0  , contact  displacement  D, 
contact  radius  tb  , and  current  particle 
radius  rp. 


tent,  and  contact  radius  must  be  known  to  use  this 
result.  It  is  unlikely  that  the  contact  radius  can  be 
determined  directly  at  any  time  for  any  given  sample, 
although  most  of  the  other  parameters  can  be  meas- 
ured. Accordingly,  it  is  necessary  either  to  know  the 
complete  history  of  a sample  of  wet  snow  or  to  re- 
late the  contact  size  to  density  through  assumptions 
about  the  geometry  of  particle  arrangements.  The 
latter  is  done  in  the  next  section,  Particle  Geometry. 

Assuming  a constant  thickness  for  the  water  film  at 
the  contact,  heat  and  water  flow  must  balance.  At  nor- 
mal rates  of  densification,  the  thickness  of  the  water 
film  at  the  contact  is  apparently  many  times  greater 
than  an  adsorbed  layer  of  structured  water.  Therefore, 
the  flux  of  water  through  the  Rim  q is  described  by 


I2p  3 r 


(6) 


where  d=  thickness  of  the  water  film 
p = pressure  in  the  liquid  phase 
p - viscosity  of  water 
r-  the  radial  coordinate. 


By  continuity,  the  flux  of  water  increases  with  radius 
along  the  contact  as 


q^^L. 

p6  dt  d 


(7) 


Hence,  for  flat  contacts,  the  pressure  distribution  along 
the  water  film, 


'"'•-Sr?'’ 

a pt  at 


(8) 


is  very  sensitive  to  the  thickness  d of  the  film.  For 
small  increases  in  the  thickness  of  the  film,  a large 
change  in  pressure  and  flux  occurs.  Therefore,  the  film 
size  must  be  relatively  stable  for  large  fluctuations  in 
heat  and  water  flow.  Unfortunately  q and  d cannot  be 
identified  independently  without  making  some  assump- 
tion about  the  pressure  field. 

Although  that  is  unnecessary  here,  it  is  important 
to  note  that  pressure  cannot  be  uniform  along  the 
inter-particle  contact.  A pressure  drop  toward  the 
outside  of  the  contact  is  necessary  for  water  flow  out 
of  the  contact  and  into  the  pore  space.  The  parabol- 
ic distribution  given  by  eq  8 shows  a fairly  uniform 
liquid  pressure  over  most  of  the  contact,  with  a rapid 
pressure  drop  near  the  edge  of  the  contact.  Thus 
the  water  film  is  capable  of  supporting  the  force  be- 
tween the  particles  with  a maximum  pressure  p0  only 
slightly  higher  than  the  average  contact  stress  Fn '* 
r ~b2.  It  is  now  apparent  that  the  tension  of  the  pore 
liquid  affects  both  the  phase  equilibrium  temperature 
of  the  pore  space  and  the  liquid  pressure  distribution 
in  the  contact.  Tire  first  effect  is  manifested  by  a 
change  in  temperature  distribution  and  rate  of  densi- 
fication,  while  the  second  influences  only  the  thick- 
ness of  the  liquid  layer  at  the  contact. 


PARTICLE  GEOMETRY 

The  rate  at  which  the  particle  centers  move  closer 
together  by  contact  melting  can  be  calculated  with 
eq  5,  but  it  is  necessary  to  specify  the  packing  geo- 
metry to  construct  a model  of  wet  snow  densification 
(see  cover  of  this  report).  Visscher  and  Bolstcrli  (1972) 
found  an  average  of  6.4  contacts  per  sphere  in  random 
packings  of  uniform  spheres,  which  is  closely  approxi- 
mated by  geometrically  regular,  cubic  packing.  A 
cubic  model  considerably  simplifies  the  complicated 
geometry  of  a random  packing,  but  does  not  adequate- 
ly represent  densification  at  porosities  greater  than 
0.476.  Higher  porosities,  which  would  have  a lesser 
number  of  contacts,  may  be  represented  by  allowing 
the  number  of  contacts  to  increase  to  six  as  the  0.476 


porosity  is  approached;  however,  some  experimental 
testing  of  this  approach  must  be  made  prior  to  its 
adoption. 

Starting  with  a cubic  packing  of  equal  spheres, 
densification  proceeds  by  melting  at  inter-particle 
contacts  and  refreezing  in  the  pore  spaces.  The  mass  of 
each  particle  is  constant,  and  we  assume  that  each  par- 
ticle remains  spherical,  except  for  its  six  contact  surfaces. 
At  any  time  during  densification,  the  depth  of  the  melted 
cap  D (see  Fig.  2)  is 

D=ro-(rl-rlt  (9) 


and  the  dry  bulk  density  p of  the  sample  is 

'Hrf)’-  (.0) 

The  intersection  of  the  spheres  causes  some  error  in 
these  equations  at  higher  densities,  but  this  error  is  less 
than  1%.  The  increase  in  particle  radius  associated  with 
the  refreezing  in  the  pore  space  is  14%  of  the  cubic 
model’s  initial  radius,  and  it  must  be  accounted  for  in 
the  model. 

Since  mass  is  conserved  for  each  particle,  the  parti- 
cle and  contact  radii  are  related  to  the  initial  radius  by 

'o  = ^p-SO-p-'b)^  -4-5'b('-p-'b),/j  • 0 1) 

In  principle  it  is  now  possible  to  solve  for  density  as  a 
function  of  time  for  any  given  set  of  conditions.  How- 
ever, algebraic  problems  make  this  approach  imprac- 
tical, and  an  iterative  calculation  of  the  deformation 
starting  from  the  cubic  packing  of  equal  spheres  is 
used  instead.  As  an  example  of  the  geometrical  rela- 
tionships and  as  an  aid  in  the  iterative  calculation,  D 
and  rb  are  shown  as  functions  of  rp  in  Figure  3. 

Results  derived  from  the  model 

As  an  example  of  the  application  of  these  concepts, 
the  hydrostatic  compression  of  saturated  snow  at  a 
hydrostatic  stress  of  1 04  N m'J  * is  calculated.  Cal- 
culated values  at  steps  in  the  iterative  procedure  are 
gven  in  Table  1 and  tire  results  are  shown  on  Figure 
4.  The  depth  of  the  melt  cap  D increases  rapidly  at 
first,  because  of  the  high  stress  concentration  at  the 
points  of  the  spherical  particles,  but  the  rate  of  growth 

* lNm"1  = IPa. 
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Figure  3.  Depth  oj  melted  cap  D and  bond  radius  rb  shown  as  a func- 
tion of  particle  radius  rp  for  a particle  with  an  initial  radius  r0  of  I mm. 


Time  (s) 

Figure  4.  The  ice  density  and  depth  of  the  melted  cap  D shown  as  fine- 
tions  of  time  for  a saturated  sample  compressed  hydrostatically  with  a 
stress  of  10*  Mm'1. 


of  the  cap  decreases  as  the  cap  volume  increases. 
Likewise,  the  rate  of  density  increase  is  a steadily  de- 
creasing function  of  time  which  approaches  the  den- 
sity of  solid  ice  asymptotically.  Qualitatively  these 
results  are  very  similar  to  those  of  Tusima  ( 1 973), 
although  his  tests  were  done  with  repeated  loadings 
at  a stress  of  3.25x  10s  N m"J  and  showed  a much 
more  rapid  response  (see  Fig.  5). 


The  rate  of  density  increase  is  a complicated  func- 
tion of  stress  because  of  the  complex  relationships 
among  the  thermodynamic  and  geometrical  variables 
which  control  the  rate  of  density  increase.  However,  as 
shown  in  Table  II,  the  rate  of  growth  of  the  melt  cap 
increasts  approximately  linearly  with  the  increase  in 
stress.  Therefore  the  rate  of  increase  in  density. 


; 


s 


' 


5 


0 667  1333  2000  2666  3333  seconds 


Figure  5.  Snow  density  shown  as  a function  of  time  for  tests  at  three  liquid 
water  contents.  The  samples  were  repeatedly  loaded  (36  loads  per  minute ) 
with  a stress  of  3.25 x 10*  N m'1.  These  curves  were  drawn  from  data  in 
Tusima  (1973). 


Table  I.  The  hydrostatic  compression  of 
saturated  snow  at  104  N m'*. 


Table  II.  The  hydrostatic  compression  of 
saturated  snow  at  two  stresses. 


rp 

rb 

D 

dD/dt 

P 

'b 

dD/dt  (7(7*  Nm 2 ) 

dD/dt  (It?  Nm'2) 

(mm) 

(mm) 

(mm) 

(mm  s l) 

(kg  m'2) 

(Nm ■*) 

(mm) 

(mm  s' ' ) 

(mm  s‘ 1 ) 

1.0 

0. 

0 

0 

480 

0 

1.0 

OO 

OO 

1.003 

0.295 

0.0413 

2,078X  10'’ 

44  5 

112,000 

1.003 

2.078XI0'7 

2.I00X  10  ‘ 

1.01 

0.4 

0.072 

I.C60X  10'’ 

602 

68,100 

1.01 

1,060X1  O'7 

1.070X  10'5 

1.022 

0.487 

0.101 

618X10'7 

662 

44,800 

1.022 

61 8X  10‘7 

622X10'5 

1.038 

0.589 

0.125 

383  XI  O'7 

718 

31,100 

1.038 

383X  10-’ 

391X  10'5 

1.062 

0.634 

0.149 

262X10'7 

780 

22,800 

1.062 

262X10'7 

268X10'5 

1.105 

0.731 

0.172 

168XI0'7 

845 

16,050 

1.105 

168XI0-7 

173X10'1 

dp  _ 3Pq'o  dD 
dt  (r0-Df  dt 
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is  approximately  linearly  related  to  stress  when  the 
load  is  first  applied,  but  is  a nonlinear  function  of 
stress  at  larger  values  of  time  because  the  denominator 
of  eq  12  decreases  with  increasing  deformation. 

The  decrease  in  phase  equilibrium  temperature  as- 
sociated with  increasing  capillary  pressure  reduces  the 
heat  flow  from  the  particle  surface  to  the  inter-par- 
ticle bond,  as  described  previously.  Thus,  at  lower 
liquid  saturations,  the  thermodynamic  deformation 
of  wet  snow  proceeds  more  slowly,  as  shown  by 


Tusima’s  (1973)  results  in  Figure  5.  Tusima’s  data 
suggest  that  the  maximum  density  approaches  a lower 
limit  at  lower  water  saturations,  an  expected  result  of 
lowering  the  phase  equilibrium  temperature  on  the 
particle’s  curved  surface.  As  a further  example  of 
the  effect  of  capillary  pressure  on  the  rate  of  defor- 
mation of  wet  snow,  the  calculated  densities  are 
shown  as  functions  of  time  at  a hydrostatic  stress 
of  104  N m"J  and  capillary  pressures  of0,0.1,0.5 
and  1 m of  water  (0,  987,  4935  and  9870  N m‘J ) 
in  Figure  6 and  Table  III.  Although  the  stress  levels 
were  higher  in  Tusima’s  experiments  where  the 
stress  was  applied  intermittently,  Tusima’s  results 
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Figure  6.  Calculated  density  shown  as  a function  of  time  for  four  values  of  capillary  pressure 
at  a hydrostatic  stress  of  H t N m3and  an  initial  particle  size  of  1 mm.  At  capillary  pressures 
of  0.5  and  1 mH20,  the  limiting  value  of  density  is  reached.  These  four  values  of  capillary 
pressure  have  the  same  effect  as  would  sodium  chloride  concentrations  ofO,  4.32 x Iff*. 
2.18x10'  and  4.36x  103  g-moles  per  kg  of  water  in  saturated  snow  at  the  same  stress. 


Table  III.  The  calculated  deformation  rates  ate  shown  for  a hydrostatic 
stress  of  104  N m * for  four  values  of  capillary  pressure. 


/*c-0 

0.1 

as 

l.OmHjO 

rp 

'b 

D 

P 

dp/dt 

dp/dt 

dp/dt 

dp/dt 

(mm) 

(mm) 

(mm) 

(Mgm3) 

(Mum'  s') 

(Mgm'J  s' ) 

(Mg  m'*  s’1 ) 

(Mg  m '*  j'1; 

1 

0 

0 

0.480 

0© 

oo 

OO 

OO 

1 + 

0.164 

0.136 

0.5046 

0.00237 

0.00231 

0.00210 

_ 

1.001 

0.221 

0.0236 

0.516 

9.64X10'* 

9.14X10'* 

7.46X  10"* 

5.36X  10"* 

1.002 

0.265 

0.0338 

0.532 

S.64X  10'* 

5.21X  10* 

3.75X  10"* 

1.93X10* 

I.OOS 

0.336 

0.0527 

0.565 

2.89X  10'* 

2.53X10* 

1.28X10-* 



1.01 

0.400 

0.0724 

0.602 

I.78X  10'* 

I.45X  10'* 

3.08X10'* 

1.02 

0.475 

0.0975 

0.653 

I.I1XI0'* 

8.00X10'* 



1.04 

0.566 

0.128 

0.724 

6.98X10'* 

3.99X10'* 

__ 



1.0S 

0.600 

0.138 

0.750 

5.99X10'* 

3.03X10'* 

_ 



1.075 

0.677 

0.157 

0.802 

4.52X10'* 

1.62X10'* 

_ 



1.10 

0.722 

0.170 

0.839 

3.66X10'* 

8.I9X  10'* 

_ 



1.14 

0.794 

0.182 

0.877 

2.80X10* 

7.96X  10'7 

_ 



1.146 

0.804 

0.183 

0.881 

2.71X10'* 

2.24X10'* 

_ 

shown  in  Figure  5 are  qualitatively  very  similar  to  the 
theoretical  results  shown  in  Figure  6.  The  rate  of 
density  increase  is  a decreasing  function  of  time,  and 
at  lower  water  saturations,  corresponding  to  larger 
capillary  pressures,  the  maximum  density  asymp- 
totically approaches  a lower  limit.  Once  this  limit  is 
reached,  further  compression  is  only  possible  by  more 
efficient  particle  packings  or  intracrystalline  deforma- 
tion. The  former  is  not  an  effective  compressive  mech- 
anism at  these  higher  densities,  while  the  latter  occurs 
much  more  slowly  than  the  thermodynamic  metamor- 
phism described  here. 
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The  introduction  of  ionic  impurities  into  the  liquid 
phase  has  a dramatic  effect  on  the  mechanical  proper- 
ties of  wet  snow  because  the  phase  equilibrium  temp- 
erature of  the  mixture  is  reduced.  When  sufficient  im- 
purities are  present,  the  mixture  temperature  can  be 
lowered  beneath  the  inter-particle  bond  temperature, 
in  which  case  the  bonds  freeze  and  a significant  inter- 
particle strength  develops.  Even  if  the  bonds  do  not 
freeze,  the  heat  flow  away  from  the  bonds  is  reduced 
and  the  rate  of  densification  and  the  maximum  density 
are  reduced  accordingly.  Thus  ionic  impurities  in  dilute 


i 


Figure  7.  Maximum  radius  of  the  inter-particle 
bond  as  a j unction  of  sodium  chloride  concen- 
tration, shown  for  an  initial  particle  size  of  1 
mm  and  a bulk  stress  of  2320  N m'5  in  satu- 
rated snow.  The  bond  size  stops  increasing 
when  the  pressure  melting  temperature  of  the 
bond  surface  equals  the  impurity  temperature 
depression  of  the  curved  particle  surface. 

concentrations  have  the  same  effect  on  mechanical 
properties  as  capillary  pressure.  To  illustrate  this  ef- 
fect, the  calculated  curves  in  Figure  6 could  be  des- 
cribed as  the  effects  of  0,4.32x  10’4,  2.18x  10‘3  and 
4.36x  10‘3  g-moles  of  sodium  chloride  per  kg  of  water 
in  saturated  snow,  instead  of  capillary  pressures  of  0, 
987, 4935  and  9870  N m‘2,  respectively.  These  levels 
of  ionic  impurities,  however,  are  not  generally  found 
in  natural  snow  covers  away  from  contaminants  such 
as  road  salt,  while  these  capillary  pressures  are  com- 
mon. 

The  effect  of  impurities  is  clearly  demonstrated  in 
Figure  7.  Like  capillary  pressure,  ionic  impurities 
limit  the  thermodynamic  densification  of  wet  snow  by 
lowering  the  phase  equilibrium  temperature  at  the 
curved  particle  surface.  For  any  given  stress  and  initial 
particle  size,  the  inter-particle  bonds  develop  to  a cer- 
tain size  before  thermal  equilibrium  is  reached.  The 
limiting  inter-particle  bond  radius  rb  is  shown  in  Fig- 
ure 7 as  a function  of  sodium  chloride  concentration. 
The  rapid  thermodynamic  densification  of  wet  snow 
is  limited  by  the  maximum  size  of  the  bond  at  thermal 
equilibrium. 


Nonhydrostatic  stress 

The  thermodynamic  deformation  of  wet  snow  is 
most  easily  described  when  a hydrostatic  bulk  stress 
causes  isotropic  densification.  In  this  case,  all  the  con- 
tacts in  the  cubic  model  experience  equal  stresses, 
temperatures,  and  rates  of  melting.  Nonhydrostatic 
states  of  stress  are  more  difficult  to  interpret,  but  some 
important  conclusions  can  be  drawn.  With  confined 
compression,  for  example,  the  lateral  stresses  are  re- 
duced by  a factor  determined  from  Poisson’s  ratio,  and 
while  the  largest  deformation  will  occur  along  the  axis 
of  compression,  the  lateral  stresses  will  give  some  lat- 
eral deformation.  In  a model  where  a constant  sample 
area  is  maintained,  a rearrangement  of  the  particles 
must  occur  to  maintain  a constant  cross-sectional  area. 
The  cubic  model  can  account  for  this  rearrangement 
only  by  allowing  the  discrete  motion  of  particles  to 
ensure  a complete  packing  of  particles  in  the  sample 
area.  Thus  the  deformation  along  the  axis  of  compres- 
sion would  be  the  sum  of  contact  melting  along  that 
axis  and  particle  displacement  from  contact  melting 
in  the  lateral  direction.  Other  states  of  stress  can  be 
treated  by  similar  techniques. 

If  one  of  the  principal  stresses  is  tensile,  the  re- 
sponse of  wet  snow  can  vary  from  immediate  failure 
at  melting  bonds  to  considerable  strength  at  frozen 
bonds.  As  described  in  the  Introduction , the  degree 
of  freezing  at  the  bonds  is  determined  by  parameters 
such  as  impurity  content,  water  saturation,  contact 
size,  and  history  of  the  sample.  If  the  bond  survives 
the  initial  tensile  loading,  the  equilibrium  temperature 
will  be  raised  by  the  tensile  stress  at  the  contact,  and 
freezing  will  strengthen  tire  contact. 

CONCLUSIONS 

The  thermodynamics  of  wet  snow  controls  its  rapid 
deformation  in  a manner  unlike  the  deformation  of 
dry  snow.  Because  of  the  presence  of  liquid  water,  wet 
snow  can  deform  rapidly  by  melting  at  the  inter-particle 
contacts  and  refreezing  at  the  curved  ice  surfaces.  Tire 
rate  of  deformation  is  controlled  by  the  parameters 
which  determine  the  temperature  difference  between 
the  inter-particle  contacts  and  the  curved  ice  surface.  The 
temperature  of  the  contact  is  determined  strictly  by 
the  force  between  the  particles  and  the  size  of  the  con- 
tact. The  temperature  of  the  curved  surface  is  determined 
by  the  particle  size,  ionic  impurity  content  of  the  liquid 


phase,  and  air  bubble  size  or  capillary  pressure  (de- 
pending on  the  saturation  regime).  As  the  inter-par- 
ticle bond  size  increases,  the  inter-particle  stress  de- 
creases and  the  bond  temperature  increases  until  the 
temperature  differences  around  a particle  are  equal- 
ized. At  that  time  rapid  thermodynamic  deformation 
ceases,  and  further  deformation  must  proceed  by 
more  efficient  packing  arrangements  and/or  the  nor- 
mal processes  of  crystalline  deformation. 

Using  a simple  geometric  packing  arrangement  and 
the  principles  of  thermodynamics,  a simple  model  of  wet 
snow  deformation  has  been  constructed.  The  results  of 
this  model  are  qualitatively  similar  to  the  experimental 
resul  ts  of  Tusima  (1 973  ),  including  a lower  value  of 
maximum  density  at  larger  capillary  pressures.  The 
effect  of  ionic  impurities  on  the  maximum  bond  size 
has  also  been  calculated  and  found  to  be  the  same  as 
increasing  the  capillary  pressure.  The  hardening  ef- 
fects of  increasing  capillary  pressure  and/or  ionic  im- 
purity content  have  practical  implications  for  strength- 
ening snow  surfaces  as  an  aid  to  oversnow  mobility. 
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